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I. 



INTRODUCTION 



The complexity involved in designing a control system 
cf reasonably low order for dynamic systems cf substantially 
higher order is a subject about which much has been written. 
The list cf references at the end of this study include only 
a small portion of the material available concerning this 
subject. A new development has been introduced by Maurice 
F. Hutton and Bernard Friedland, [Ref 1], which provides a 
systematic and analytical approach to obtaining 
reduced-order approximations from higher-order transfer 
furcticns. This method is referred to as the "South 
Approximation Method," and is based on an expansion that 
uses the South table of the original higher-order transfer 
function. The desire for low order models to simulate 
higher-order practical systems, such as electrical power 
plants, chemical processes, aircraft model designs and 
electronic circuitry is obvious when dealing with such 
complex systems. The "Routh Approximation Method", is thus, 
an extremely powerful tool to utilize in dealing with these 
types of prchlems. 

In comparison to other available methods presently 
used, the South method has many advantages. The widely used 
"Dominant Pole", approximation method , which is based on 
approximating a system by utilizing the poles, or 
characteristic roots nearest the imaginary axis, has the 
disadvantage that the roots of the characteristic equation 
must be found. For a very- high-order system, this is not a 
trivial task. The "Pad'e" approximation is based on setting 
the numerator and denominator orders tc a desired value, the 
coefficients are then chosen so that the Taylor series 
expansions cf the approximant and the original transfer 
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function agree in as many terms as possible. This method 
produces accurate results, however, it is limited in 
application to single-input, single-output system analysis. 
In addition, an unstable approximant may be obtained frcm a 
stable system since the approximant ' s poles depend on both 
the original equation's numerator and denominator. 

The "South Approximation Method," preserves stability 
if the original transfer function is stable. It provides an 
efficient means of obtaining lower-order approximants for 
multiple inputs or octputs and is very adaptable to computer 
programming. 



This thesis is concerned with the development of a 
computer program which utilizes the "Routh Approximation 
Method," to obtain lower-order transfer functions from 
higher-order functions and compares the resulting eguations 
to the original eguations by graphically displaying the 
response of each to various input excitations. The program 
was designed to provide output data which is useful in 
determining which degree of approximant is best suited to 
simulate the original higher-order equation. The orders 
available are the first through the fourth order reduced 
transfer functions. The program is capable of reducing 
transfer functions up to, and including the tenth order and, 
without loss of generality, may be extended to handle an 
almost unlimited order. 



For illustrative purposes, a seventh order system was 
used, as an example, to demonstrate the computer program's 
capabilities. This particular system was chosen merely to 
demonstrate the simplicity involved in utilizing the program 
and has no physical relationship to any real system. 
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II. NATURE CF THE PROBLEM 



In approximating higher-crder systems by lower-order 
models, a linear system is desirable since the linear 
characteristic eguations are less complicated than the 
non-linear eguaticns. Therefore, this study is primarily 
concerned with linear, time invariant characteristic 
eguations, since the objective in developing the computer 
pregram was tc illustrate a simple analytical approach to 
obtaining lewer erder char acteristic eguations from higher 
order systems. 

It is obvious that a reduced-order model cannot 
characterize a given system as accurately as a higher-order 
model. The validity of the lower-order model is based upon 
its degree of success in approximating the higher order 
system in representing the characteristics of primary 
interest . 

Interpreting the solution of a higher-order system 
often results in computational difficulties which are 
reduced by appropriate selection of a reduced-order 
ap p reximatien . 

Ideally, a reduced-order model would approximate the 
higher-orcer system in both lew and high frequency ranges. 
In doing sc, seme accuracy is lost in compensating fer the 
different responses of the system to variations in 
frequency . 

The lew frequency model more closely approximates the 
higher-order system than a model composed cf both lew and 
high erder frequency characteristics. The procedure 
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utilized in this study places emphasis on the low frequency 
model . 

For unstable systems, the program described herin is 
still valid, however, the original higher-order system must 
first be modified by shifting the imaginary axis prior to 
the approximation. Thus for an unstable transfer 
function , 8 (s) , the equation must be changed tc H (s+a) , where 
a is chcser sufficiently large so that H(s) is 
asymptotically stable. This procedure is described in 
detail in [Eef 1], pp 332. 
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III. COMPOTES PROGRAM CRITERIA 



A. GENE 5 AL 

The program, referred to as SOUTH'!, was written with 
the following criteria: 

1 . Mini mum u til i z atio n of co m pu ter time 

R00TH1 consists of less than 150K of storage and takes 
less than 2 minutes of computer time. This net only 
conserves efficiency, but also provides the user with the 
desired data at a minimum cost. 

2 . Ease cf use 

Input data required consists of eleven data cards for 
maximum utilization of the program. Emphasis is placed on 
the ease cf using the program to obtain the desired results 
with minimum time expended on computer programing. 

Usef u 11 n ess of ou tpu t 

a. The first, second, third and fourth order 
approximants to the original equation are printed in 
transfer function format. Both numerator and denominator 
coefficients are printed in ascending powers of S. 

b. The roots of the original and reduced equations 
are provided to enable the user to study the response cf the 
systems in the frequency domain. 
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c. Choice of variables for print out in table form 
is available with up to eight variables maximum for any one 
cf three possible runs. 



d. Choice of variables for graphical output with up 
to four curves may be plotted separately or all cn one 
graph . 



consisting 
available, 
compared 
plotted tc 
inputs may 



Graphical output response to input excitations 
of a Step, Ramp, or Sinusoidal input are 
The graphs display the original output response 
to the lower order response and the error is 
display the differences in response. Multiple 
be used. 



B. SUBROUTINES 



Two subroutines are utilized in R0UTH1. The rocts of 
the original and reduced equations are determined by 
subroutine PRQD , »hicn was taken from [Ref 3], and the 
tables and graphs are determined by subroutine R2DUCT1 which 
is a modification of INTEGl , [Ref 3]. 
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IV- FOURTH ORDER EXAMPLE 



A. THREE STEP PROCEDURE 

To illustrate the "Routh Approximation Method", A 
fourth order example was chosen for simplicity. 

The method, described in detail in [Ref 1], consists of 

three basic steps. The following transfer function is used 
to illustrate the procedure : 

H (s ) = * 

20 + 32S + 24^ + 8S 5 + S 4 

The first step is to compute what are termed Alpha and 
Beta coefficients from the Routh Table shown on page 19. 
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Step two in the procedure is to obtain what are termed 



the Routh convergents , which are based on the following : 

Letting A^s) and B k (s) denote the denominator and the 
numerator, respectively, of the Routh convergent, i.e., 

A^(s) = Alpha^S + 1 
(s ) = Beta^ 

2 

A 2 (s) = Alpha^Alpha 2 S + Alpha 2 S + 1 
B 2 (s) = Alpha 2 Beta^S + Beta 2 



The general expression from [Ref 1] is the following : 
A k (s) = Alpha k (s) A k-1 (s) + A k _ 2 (s) 

B k (s) = Alpha k SB k _ 1 (s) + B k _ 2 (s) + Beta k k=l,2,3... 

with A_^(s)=0.0 B_^(s)=0.0 

A q ( s ) =1.0 B Q (s) ”°*° 

The Routh convergents for the fourth order example are 
the following : 

Ri (s)— 0.0312S 
0 . 625S + 1 

R 2 (s)= 0.056316 ' 

1.05263S 2 + 1.6842S + 1.0 

R 3 (s)= 0 . 158833S 2 — 0.0083 

3.166S 3 + 5.066S 2 + 3.633S + 1 

R 4 (s)= sf 

20S 4 + 32S 3 + 24S 2 + 8S + 1.0 
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The third, and final step in the procedure is to apply 
what is termed a reciprocal transformation, defined by 

H k (s) = s xR k (1/s) 

which is merely a reversal of the order of the polynomial 
coefficients. Thus, for the example given, the reduced order 
approximations are given by the following : 

H, (s) = 0.03125 

S + 0.625 

H 2 (s) = 0.056316 

5 2 + 1.6842S + 1.05263 

H 3 (s) = -0.0083S 2 + 0.158333 

5 3 + 3.633S 2 + 5.066S + 3.166 

H 4 (s) = 1 

5 4 + 8S 3 +24S 2 + 32S + 20 

/ 

As expected, the fourth-order approximation is the same 
as the original equation. 

The poles of the approximants are illustrated in table 
IV. 2. The poles of the approximants approach the dominant 
poles of the original higher-order equation as the order of 
the approximant is increased. 
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Poles and Zeros of Approximant 


Order 


Poles 


Zeros 


R 2 


-0. 842 + JO .586 


0 


R 3 


-0.920 + JO. 958 
-1.792 


+ 4.36 


R 4 


-2.0 + J1.0 
-2.0 
-2.0 


0 


A 


-2.0 + J1.0 

-2.0 

-2.0 


0 




+ 



-4 




POLES OF ROUTH APPROXIMANTS 
Table IV. 2 
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V 



COMPUTES EXAMPLE 



A. SEVENTH CEDES SYSTEM 



A seventh order system with the transfer function given 
by [ Ref 2 ] / as 

,7 



384x10 



G ( s ) — — j 



S'+432S 6 +62670S 5 +3615900S 4 +75114000S 3 +553920000S 2 . 



.+1443200000S+3 84x10 



or in factored form as 



G(s) = 



384x10 



(S +2S+10) (S+10) (S+20 ) (S+80 ) (S+120) (S+200) 

was reduced to the low order approximan ts by 
utilizing ROUTH1 and the response to input excitations, 
consisting of a Step and Ramp are illustrated in figures 5.1 
through 5.5. 



The rocts of the system and its low order apprcximants 
are illustrated in table V.l. 
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Poles and Zeros of Approximant 


Order 


Poles 


Zeros 


R 2 


-2.038 + j 2 . 5 86 


0 


R 3 


-1.084 + '-j'ZTVTr 
-6.291 


+ 10.132 


R 4 


-1.0 + j 2 . 999 
-10.89 + j 2 . 30 


+ 316.58 


A 


-1.0 + j 2 . 999 
-10.0 
-200.0 
-120.0 
-80.0 

^22^2 


0 



© 




©■ 



*3 



Im 



© 




1 



Re 



-10 



© 



-8 



-6 



-4 



-2 



..-1 



---2 



Q 



POLES OF ROUTH APPROXIMANTS 
Table V.l 
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Figure 5 . 1 
STEP INPUT 
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Figure 5 . 2 
STEP INPUT 



24 



MG i’TC C’C SCO ‘(CO 'CO ?X 




Figure 5 . 3 
STEP INPUT 



25 



o 

a 




Figure 5 . 4 
RAMP INPUT 
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Figure 5 . 5 
RAMP INPUT 
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Figure 5 . 6 
RAMP INPUT 
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VI. COMPARISON TO OTHER METHODS 



A. DOMINANT POLE METHOD 

The Dcainant Pole Approximation method is oased upon 
utilizing the poles closest to the imaginary axis. The 
equation must be factored to obtain the characteristic 
roots. Fcr the given seventh-order system, the lower- order 
equations are given by this method as the following: 
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H (s ) =- 



S +2S+10 



100 



H(s) = 



S 3 +14S 2 +30S+100 



2000 



H ( S ) = 



S 4 +34S 3 +310S 2 +700S+2000 



Graphical plots of the Dominant Pole reduced 
equations, in response to Step inputs, are illustrated in 
figures 6.1 through 6.3 in comparison to the original 
equation and figures 6.4 through 6.6 illustrate the 
comparison to the South equations. The error between the 
systems is also plotted. Table VI. 1 gives the analytical 
data for performance measure comparisons. 
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B 



ITERATIVE OPTIMIZATION METHOD 



The following equations, representing the seventh-order 
example, were taken from [Ref 2]. This method makes use of 
an iterative ninimization technique to locate the best pole 
and zero locations for the lower-order models. 

7.203856 

H (s ) = — ~ 

S+1.98616S+7. 203856 



H (s 



52.2861 

S 3 +7.1383S 2 +19.5015S+52.2861 



1470.1403 

H (s) =— 4 3 2 

S +28.5204S J +209.3842S +552.5241S+1470.1403 



Plots of these equations v 
equation and the Routh equations are 
6.7 through 6.12 in response to Step 
provides numerical data for performanc 



ersus the Original 
illustrated in figures 
inputs. Table VI . 1 
e measure comparisons. 
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C. PRESENTATION OF DATA 

Comparisons were made between the "Routh Approximation 
Method", and that of the Dominant Pole and Iterative 
Optimization methods previously described. 

The basis for comparison consists of the following 
criteria: 

1. Peak Overshoot M =Y^ 

pt max' ss 

2. Delay Time T =time for Y(t) to reach 0.5 Y 

J d ss 

the very first time. 

3. Rise Time T r =time for Y(t) to go from 0.1 to 

0.9 of the final value. T r =l/BW. 

4. Settling Time T =time at which Y(t)=Y 

s s s 

5. Graphical Representation in response to Step inputs. 



6 . 




t^=integration steps 



Table VI. 1 illustrates the above comparisons 
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Figure 6 .1 
STEP INPUT 
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Figure 6 . 3 
STEP INPUT 
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Figure g . 4 

STEP INPUT 
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Figure 6 .5 
STEP INPUT 
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Figure 6 . 6 
STEP INPUT 
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Figure 6 . 7 
STEP INPUT 
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Figure 5 . 8 
STEP INPUT 
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Figure 6.9 
STEP INPUT 
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Figure 6 10 
STEP INPUT 
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Figure 6 . 11 
STEP INPUT 
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Figure 6 . 12 
STEP INPUT 
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TABLE VI. 1 



Table A = Original Seventh Order 

Symbols R _ R OU -|-h Approximation 

D = Dominant Pole Approximation 

O = Optimum Minimization Method 



PERFORMANCE MEASURE COMPARISONS 
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VII. CONCLUSION 



The graphical data and tabulated analysis of this 
thesis indicate that the South Approximation Method is a 
valuable formulation technique which produces very 
satisfactory results in acquiring approximations to 
higher-order systems with minimum cost. 

In comparison to the other methods of analysis, the 
Rcuth Approximation Method offers a guick and easy 
analytical approach to obtaining lew order models. Unlike 
the Fad'e Approximation, the South method ensures stability 
of the lewer-erder models, if the original higher-order 
system is stable. The original system need net be factored, 
as in the Dcninant Pole Approximation method. 

The computer program, H0UTH1, utilized to acquire the 
reduced order equations, the roots of the eguations and the 
graphical plots and numerical tables, takes considerably 
less time than the minimization technique utilized in 
reference 2. However, the minimization technique operates 
efficiently without prior knowledge of the system’s transfer 
function. In the South method the availability of the 
transfer function is a necessary requirement. 

In comparing the low order equations of the various 
methods discussed, the Routh method has proven to be a valid 
and efficient solution to the problem of obtaining good low 
order approximants to complex higher-order systems. 
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